New relativistic dissipative fluid dynamics from kinetic theory 
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Starting with the relativistic Boltzmann equation where the collision term is generalized to include 
nonlocal effects via gradients of the phase-space distribution function, and using Grad's 14-moment 
approximation for the distribution function, we derive equations for the relativistic dissipative fluid 
dynamics. We compare them with the corresponding equations obtained in the standard Israel- 
Stewart and related approaches. Our method generates all the second-order terms that are allowed 
by symmetry, some of which have been missed by the traditional approaches based on the 14- 
moment approximation, and the coefficients of other terms are altered. The first-order or Navier- 
Stokes equations too get modified. Significance of these findings is demonstrated in the framework 
of one-dimensional scaling expansion of the matter formed in relativistic heavy-ion collisions. 



I. INTRODUCTION 

Relativistic fluid dynamics finds applications in cos- 
mology, astrophysics, and the physics of high-energy 
heavy-ion collisions. In cosmology and certain areas of 
astrophysics, one needs a fluid dynamics formulation con- 
sistent with the General Theory of Relativity [l[ . On the 
other hand, a formulation based on the Special Theory of 
Relativity is quite adequate to treat the evolution of the 
strongly interacting matter formed in high-energy heavy- 
ion collisions when it is close to a local thermodynamic 
equilibrium. The correct formulation of the relativistic 
dissipative fluid dynamics is far from settled and is cur- 
rently under intense investigation @-H|. 

In applications of fluid dynamics it is natural to first 
employ the zeroth order (in gradients of the hydrody- 
namic four-velocity, for example) or ideal fluid dynam- 
ics. However, as all fluids are dissipative in nature due to 
the uncertainty principle (To| , the ideal fluid results serve 
only as a benchmark when dissipative effects become im- 
portant. The first-order dissipative fluid dynamics or 
the relativistic Navier- Stokes (NS) theory [Tl| involves 
parabolic differential equations and suffers from acausal- 
ity and instability. The second-order Israel-Stewart (IS) 
theory [l2| , with its hyperbolic equations restores causal- 
ity but may not guarantee stability (l3j . 

The second-order viscous hydrodynamics has been 
quite successful in explaining the spectra and azimuthal 
anisotropy of particles produced in heavy-ion collisions at 
the Relativistic Heavy Ion Collider (RHIC) [H, JT1 and 
recently at the Large Hadron Collider (LHC) 16. [l7| 



However, IS theory can lead to unphysical effects such as 
reheating of the expanding medium [181 ] and to a negative 
pressure [i~9j at large viscosity indicating its breakdown. 
Furthermore, from comparison to the transport theory 
it was demonstrated [f| [l3[ that IS approach becomes 
marginal when the shear viscosity to entropy density ra- 
tio rj/s > 1.5/(47r). With this motivation, the dissipative 
hydrodynamic equations were extended [7| to third order, 
which led to an improved agreement with the kinetic the- 
ory even for moderately large values of rj/s. 

It is well known that the approach based on the gen- 
eralized second law of thermodynamics fails to capture 



all the terms in the evolution equations of the dissipative 
quantities when compared with similar equations derived 
from transport theory [2j]. It was pointed out that using 
directly the definitions of the dissipative currents, instead 
of the second moment of the Boltzmann equation as in 
IS theory, one obtains identical equations of motion but 
with different coefficients [8[ . Recently, it has been shown 
9f| that a generalization of Grad's 14-moment method 
20| results in additional terms in the dissipative equa- 
tions. 

It is important to note that all formulations that em- 
ploy the Boltzmann equation make a strict assumption 
of a local collision term in the configuration space [8j, [l2j . 
In other words, within an infinitesimal fluid element con- 
taining a large number of particles and extending over 
many interparticle spacings llf, the different collisions 



that increase or decrease the number of particles with a 
given momentum p are all assumed to occur at the same 
point x^ . This makes the collision integral a purely local 
functional of the single-particle phase-space distribution 
function f(x,p) independent of the derivatives d tl f(x,p). 
In kinetic theory, f(x,p) is assumed to vary slowly over 
space-time, i.e., it changes negligibly over the range of 
interparticle interaction [211 ]. However, its variation over 
the fluid element may not be insignificant; see Fig. [1] 
Inclusion of the gradients of f{x,p) in the collision term 
will affect the evolution of dissipative quantities and thus 
the entire dynamics of the system. 

In this Letter, we shall provide a new formal derivation 
of the dissipative hydrodynamic equations within kinetic 
theory but using a nonlocal collision term in the Boltz- 
mann equation. We obtain new second-order terms and 
show that the coefficients of the other terms are altered. 
These modifications do have a rather strong influence on 
the evolution of the viscous medium as we shall demon- 
strate in the case of one-dimensional scaling expansion. 



II. NONLOCAL COLLISION TERM 

Our starting point is the relativistic Boltzmann equa- 
tion for the evolution of the phase-space distribution 
function, p^d^f = C[f], where the collision term C[f] 



is required to be consistent with the energy-momentum 
and current conservation. Traditionally C[f] is also as- 
sumed to be a purely local functional of f(x,p), inde- 
pendent of du f ■ This locality assumption is a powerful 
restriction 12] which we relax by including the gradients 
of f(x,p) in C[f], This necessarily leads to the modified 
Boltzmann equation 

P^J = C m [f] = C[f] + dniAPf) + d^(B^f) + ..., 

(1) 

where A^ 1 and B^ v depend on the type of the collisions 
(2 o 2, 2o 3, ...)■ 
For instance, for 2 o 2 elastic collisions, 

C[f] = ^ f dp'dk dk' W pp >^ kk , 

X(fkfk' fpfp' - fpfp' fkfk'), (2) 

where W pp >^kk' is the collisional transition rate, f p = 
f(x,p) and f p = l- rf(x,p) with r = 1,-1,0 
for Fermi, Bose, and Boltzmann gas, and dp = 
gdp/[(2n) 3 y/p 2 + to 2 ], g and m being the degeneracy 
factor and particle rest mass. The first and second terms 
in Eq. @ refer to the processes kk' — > pp' and pp' — > kk', 
respectively. These processes are traditionally assumed 
to occur at the same space-time point x^ with an under- 
lying assumption that f(x,p) is constant not only over 
the range of interparticle interaction but also over the en- 
tire infinitesimal fluid element of size dR, which is large 
compared to the average interparticle separation [TTj] ; see 
Fig. [U Equation (1) together with this crucial assump- 
tion has been used to derive the standard second-order 
dissipative hydrodynamic equations [1, [H, [Hj] . We, how- 
ever, emphasize that the space-time points at which the 
above two kinds of processes occur should be separated 
by an interval < dR within the volume d A R. It 
may be noted that the large number of particles within 
d A R collide among themselves with various separations 
£ M . Further, £ M is independent of the arbitrary point x^ 1 
at which the Boltzmann equation is considered, and is 
a function of {p' , k, k'). Of course, the points (x M — £ M ) 
must lie within the past light-cone of the point x M (i.e., 
£ 2 > and £° > 0) to ensure that the evolution of f{x,p) 
in Eq. ([1} docs not violate causality. With this re- 
alistic viewpoint, the second term in Eq. ([2J involves 
f(x - €,p)f(x ~ £,p')f(x - £,k)f{x - €,k'), which on 
Taylor expansion at x^ up to second order in results 
in the modified Boltzmann equation ([lj with 

^ = \ j dp'dkdk' eW w ^kwfv'fkfw, 

= ~\J dp'dk dk' eCWpp^wU.fkh'- (3) 

In general, for all collision types (2 2, 2 «-» 
3,...), the momentum dependence of the coefficients 
and B^ can be made explicit by expressing them 
in terms of the available tensors p^ and the metric 



INFINITESIMAL VOLUME ELEMENT IN A FLUID 




FIG. 1: Collisions kk' — > pp' and pp —¥ kk' occurring at 
points x^ and x^ — £ M within an infinitesimal fluid element 
of size dR, around a;'', containing a large number of particles 
represented by dots. 

g» v = diag(l, -1,-1,-1) as A^ = a(x)p» and B^ = 
bi(x)g f11 ' + b-2{x)p^p" ', in the spirit of Grad's 14-moment 
approximation. Equation ([1} forms the basis of our 
derivation of the second-order dissipative hydrodynam- 
ics. 



III. HYDRODYNAMIC EQUATIONS 

The conserved particle current and the energy- 
momentum tensor are expressed as [2l| 

N" = J dp p»f, T» v = f dp p»p v f. (4) 

The standard tensor decomposition of the above quanti- 
ties results in 

T" u = euV- (P + njA^+Tr"", (5) 

where P, n, e are respectively pressure, number density, 
energy density, and A Ml/ — ~ u^u v is the projec- 
tion operator on the three-space orthogonal to the hydro- 
dynamic four-velocity defined in the Landau frame: 
T^ u u v = eu M . For small departures from equilibrium, 
f(x,p) can be written as / = /o + Sf. The equilibrium 
distribution function is defined as /q = [exp(/3u -p — a) + 
r] _1 where the inverse temperature /3 = 1/T and a = j3p 
(/i being the chemical potential) arc defined by the equi- 
librium matching conditions n = uq and e = eq. The 
scalar product is defined as u.p = u^p^ . The dissipative 
quantities, viz., the bulk viscous pressure, the particle 
diffusion current and the shear stress tensor are 

n = -^f J dpp^sf, 

n" - A"" fdppuSf, 

tt"" = A^ J dpp a P ^6f. (6) 

Here = [A£A£ + A£A£ - (2/3)A^A^]/2 is the 
traceless symmetric projection operator. Conservation 



3 



of current, d^N^ = and energy-momentum tensor, 
Q^X^ V = o, yield the fundamental evolution equations 
for n, e and u M 



expansion scalar 9 = d- u, we finally obtain the following 
evolution equations for the dissipative fluxes defined in 
Eq. ©: 



Dn 



De + (e + P + n)a M u M - tt^ Vfrti* 



(e + P + Il)Du a - V"(P + n) 



0, 

o, 

0. (7) 



We use the standard notation A {a B f3) = (A a B^ 
A^B a )/2, D 



u^df,, and V a 



A MQ 9 M . For later use 



we introduce = A^X V and X^ = A^X a ? . 

Conservation of current and energy-momentum implies 
vanishing zeroth and hrst moments of the collision term 
C m [f] in Eq. ©, i.e., JdpC m [f}=0 = J dp p"C m [f\. 
Moreover, the arbitrariness in £ M requires that these con- 
ditions be satisfied at each order in Retaining terms 
up to second order in derivatives leads to three constraint 
equations for the coefficients (a, b\, b 2 ), namely d^a = 0, 



d 2 (6i(l)o) + d^(b 2 {p^p v )o) = 0, 
uadpdu (M2>W*)o) + u a d 2 (hnu a ) = 0, 



(8) 



where we define (• • • )o = / dp(- ■ ■ )/o- It is straightfor- 
ward to show using Eq. @ that the validity of the sec- 
ond law of thermodynamics, d^s^ > 0, enforces a further 
constraint \a\ < 1, on the collision term C m [f]. 

In order to obtain the evolution equations for the dis- 
sipative quantities, we follow the approach as described 
by Denicol-Koide-Rischke (DKR) in Ref. 0. This ap- 
proach employs directly the definitions of the dissipative 
currents in contrast to the IS derivation which uses the 
second moment of the Boltzmann equation. The comov- 
ing derivatives of the dissipative quantities can be written 
from their definitions, Eq. (|6]), as 



n 



• (IV 



A 



a0 



dp P a p?Sf, 



A"" / dp Pv Sf, 



dp p a P 5f, 



(9) 



where, X = DX. Comoving derivative of the nonequi- 
librium part of the distribution function, Sf, can be ob- 
tained by writing the Boltzmann equation (|T|) in the 
form, 

Sf = -fo-— p"V M / + — C m [f]. (10) 
u.p u.p 

To proceed further, we take recourse to Grad's 14- 
moment approximation (20j for the single-particle dis- 
tribution in orthogonal basis [8[ 

/ = fo + /o/o (A n n + X n n aP a + K-K^p?) . (11) 

The coefficients (An,A„,A„.) are functions of (n, e,/3, a). 
Using Eqs. ([§))- (jTTj) and introducing first-order shear ten- 
sor a^v = V^it„>, vorticity = (V^it,, — V v u^)/2 and 



n 



n 



T' Qn n ■ it - l' nn d ■ n - 5' nn m 



- \' Un n ■ Va + An 7r 7r M ^cr MI/ + A nil ii ■ u 



AnuOVW^ + (8 terms), 

T 



A( 



+ Kuu^uv + A^A^a; 7 " + (9 terms), (13) 



2AL r 7r^w">" + AL,n<"V>« - n,,-.'^ 



1 _</^}p 

7T7T " p 



+ Xih^ + X 2U { ^ u) b 2 + X3VW>&2. (14) 

The "8 terms" ("9 terms") involve second-order, linear 
scalar (vector) combinations of derivatives of 61 , 62 ■ All 
the terms in the above equations are inequivalent, i.e., 
none can be expressed as a combination of others via 
equations of motion 23]. All the coefficients in Eqs. (fT2")) - 
(|14[) are obtained as functions of hydrodynamic variables. 
For example, some of the transport coefficients related to 
shear are 



= a + [((u.p) 3 ) - m 2 n] , 



(15) 



where a = (1 — a). The rest of the coefficients will be 
given in [24|. 

Retaining only the first-order terms in Eqs. (IT2^) - (IT4l . 
and using DKR values of bulk viscosity £, thermal con- 
ductivity k and shear viscosity 77, we get the modified 
first-order equations for bulk pressure II = — r n /3 n # = 
— a£6>, heat current — (3' n T' n \/^a and shear stress ten- 
sor ir^ = 2t'^(3'^ u = lar^G^ = 2r}aa> w . Thus the 
nonlocal collision term modifies even the first-order dissi- 
pative equations. This constitutes one of the main results 
in the present study. 

If a, b\ and 62 are all set to zero, Eqs. (|12[) - (fT4"|) reduce 
to those obtained by DKR [8j with the same coefficients. 
Otherwise coefficients of all the terms occurring in the 
DKR equations get modified. Furthermore, our deriva- 
tion results in new terms, for instance those with coef- 
ficients Afc«, Aku, (k — n,n, 7r), which are absent in [8| 
as well as in the standard Israel- Stewart approach [ijj. 
Hence these terms have also been missed so far in the 
numerical studies of heavy-ion collisions in the hydro- 
dynamic framework [13, El HI . Indeed Eqs. (THl - fTI)) 
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FIG. 2: Time evolution of (a) temperature, shear pressure, 
inverse Reynolds number and parameters (&i, 62) normalized 
to their initial values, and (b) anisotropy parameter Pl/Pt- 
Initial values are to = 0.9 fm/c, To = 360 MeV, 77/s = 0.16, 
7ro = 477/ (3ro). Units of hi are GeV -2 . The curve labelled 
DKR is obtained by setting a = 61 = 62 = in Eqs. (|16[) and 

G3- 



contain all possible second-order terms allowed by sym- 
metry considerations 23[. This is a consequence of the 
nonlocality of the collision term C m [/] . However, we note 
that a generalization of the 14-moment approximation is 
also able to generate all these terms as shown recently in 
Ref. d. 



IV. NUMERICAL RESULTS 

To demonstrate the numerical significance of the new 
dissipative equations derived here, we consider evolu- 
tion of a massless Boltzmann gas, with equation of state 
e = 3P, at vanishing net baryon number density in the 
Bjorken model (26[. The new terms, namely ii-ii, w^w"* 1 , 

u}^ v ii v , A^d^uj 11 ' , ii^y,") and ujj^ lo^ p containing acceler- 
ation and vorticity do not contribute in this case. How- 
ever, they are expected to play an important role in the 
full 3D viscous hydrodynamics. 

In terms of the coordinates (r, x, y, rf) where r = 



Vt 2 



and 77 = tanh (z/t), the initial four- velocity 



becomes u M = (1,0,0,0). In this scenario n = = 
and the equation for it = — T 2 7^' ,,, reduces to 

7T „ dm „ 4 , 7T , dbn 

— + P*T = P*Tr-*--ifrr-r, (16) 

7V or or t dr 

where the coefficients are 

, - , h(e + P) 4 9(e + P) 

ft = r-5 , p 7r = -aP, tp = — — — , 

aprj 5 bapr] 



A = — a 
21 



5 



7/3 



arj 



(17) 



For comparison we quote the IS results [T2|: /3 ff = 
2P/3, A = 2. The coupled differential equations (J7J, 
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(|8]) and (|16j) are solved simultaneously for a variety of 
initial conditions: temperature T = 360 or 500 MeV cor- 
responding to typical RHIC and LHC energies, and shear 
pressure n = or ir = 7Tns = 4?7/(3to) corresponding to 
isotropic and anisotropic pressure configurations. Since 
the nonlocal effects embodied in the Taylor expansion ([1]) 
are not large, the initial a, b±, 62 are so constrained that 
the corrections to first-order and second-order terms re- 
main small; recall also the additional constraints \a\ < 1 
and Eq. ©. 

Figure 2(a) illustrates the evolution of these quantities 
for a choice of initial conditions. T decreases monoton- 
ically to the crossover temperature 170 MeV at time 
r ~ 10 fm/c, which is consistent with the expected life- 
time of quark-gluon plasma. Parameter a is constant 
whereas b\ and 62 vary smoothly and tend to zero at large 
times indicating reduced but still significant presence of 
nonlocal effects in the collision term at late times. This is 
also evident in Fig. 2(b) where the pressure anisotropy 
Pl/Pt = (P — v)/(P + 7r /2) shows marked deviation 
from IS, controlled mainly by a. At late times Pl/Pt is 
largely unaffected by the choice of initial values of b\, &2- 
Although the shear pressure tt vanishes rapidly indicat- 
ing approach to ideal fluid dynamics, the Pl/Pt is far 
from unity. Faster isotropization for initial a > may be 
attributed to a smaller effective shear viscosity (1 — a)r] 
in the modified NS equation, and conversely. Figure 2(b) 
also indicates the convergence of the Taylor expansion 
that led to Eq. Q. 

Figure 3 shows the evolution of Pl/Pt for isotropic 
initial pressure configuration, at various rj/s for the LHC 
energy regime. Compared to IS, DKR leads to larger 
pressure anisotropy. Further, with small initial correc- 
tions (10% to first-order and ~ 20% to the second-order 
terms) due to a, b±, 62, the nonlocal hydrodynamics 
(solid lines) exhibits appreciable deviation from the (lo- 
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cal) DKR theory. The above results clearly demonstrate 
the importance of the nonlocal effects, which should be 
incorporated in transport calculations as well. Compari- 
son of nonlocal hydrodynamics to nonlocal transport the- 
ory would be illuminating. 

In a realistic 2+1 or 3+1 D calculation, one has to 
choose the thermalization time and the freeze-out tem- 
perature together with suitable initial conditions for hy- 
drodynamic velocity, energy density, shear pressure as 
well as for the nonlocal coefficients a, b%, 62 to fit dN/dr] 
and pt spectra of hadrons, and then predict, for example, 
the anisotropic flow v n for a given rj/s. Nonlocal effects 
(especially via a) will affect the extraction of rj/s from fits 
to the measured v n . It may also be noted that although 
(local) viscous hydrodynamics explains the gross features 
of 7r~ and K~ spectra for the (0-5)% most central Pb-Pb 
collisions at ^/snn = 2.76 TeV, it strongly disagrees with 
the measured p spectrum [27| . Further the constituent 
quark number scaling violation has been observed in the 
v-i and vz data for p, at this LHC energy [28| . The above 
discrepancies may be attributed partly to the nonlocal 
effects which can have different implications for two- and 
three-particle correlations and thus affect the meson and 
baryon spectra differently. 



troducing a nonlocal generalization of the collision term 
in the Boltzmann equation. The first-order and second- 
order equations are modified: new terms occur and coef- 
ficients of others are altered. While it is well known that 
the derivation based on the generalized second law of 
thermodynamics misses some terms in the second-order 
equations, we have shown that the standard derivation 
based on kinetic theory and 14-moment approximation 
also misses other terms. The method presented here is 
able to generate all possible terms to a given order that 
are allowed by symmetry. It can also be extended to 
derive third-order hydrodynamic equations. 
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